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Abstract: In this paper, a new tensor-based subspace approach is proposed to estimate the 
direction of departure (DOD) and the direction of arrival (DOA) for bistatic multiple-input 
multiple-output (MIMO) radar in the presence of spatial colored noise. Firstly, the received 
signals can be packed into a third-order measurement tensor by exploiting the inherent 
structure of the matched filter. Then, the measurement tensor can be divided into two 
sub-tensors, and a cross-covariance tensor is formulated to eliminate the spatial colored 
noise. Finally, the signal subspace is constructed by utilizing the higher-order singular 
value decomposition (HOSVD) of the cross-covariance tensor, and the DOD and DOA can 
be obtained through the estimation of signal parameters via rotational invariance technique 
(ESPRIT) algorithm, which are paired automatically. Since the multidimensional inherent 
structure and the cross-covariance tensor technique are used, the proposed method provides 
better angle estimation performance than Chen's method, the ESPRIT algorithm and the 
multi-SVD method. Simulation results confirm the effectiveness and the advantage of the 
proposed method. 

Keywords: MIMO radar; DOD and DOA estimation; spatial colored noise; higher-order 
singular value decomposition 
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1. Introduction 

Recently, multiple-input multiple-output (MIMO) radar [1-3] has drawn increasing attention and 
has become a hot research topic in the area of radar. MIMO radar uses multiple antennas to emit 
simultaneously orthogonal waveforms instead of the coherent waveforms, which are used in the 
phased-array radar, and this waveform diversity endows MIMO radar with superior performance relative 
to phased-array radar. Based on the configuration of transmit and receive antennas, MIMO radar can be 
grouped into two classes. One is called statistical MIMO radar [2], which can solve the problem of target 
scintillation, due to the widely spaced transmit/receive antennas. The other is called colocated MIMO 
radar [3], including bistatic and monostatic MIMO radar [4,5], whose transmit antennas and receive 
antennas are close spaced. The colocated MIMO radar can obtain the virtual aperture, which is larger 
than the real aperture, so it brings a lot of advantages, such as narrower beamwidth and lower sidelobes, 
higher angular resolution and angular estimation accuracy. 

Angle estimation is an important aspect in array signal processing and MIMO radar [6-14]. In bistatic 
MIMO radar, the direction of departure (DOD) and the direction of arrival (DOA) need to be estimated 
simultaneously. In [8], a two-dimensional Capon estimator is applied to estimate DOD and DOA, which 
are paired automatically. However, it has a heavy computational burden owing to the two-dimensional 
spectrum searching. In order to alleviate the computational burden, the estimation of signal parameters 
via rotational invariance technique (ESPRIT) [9,10] is employed to DOD and DOA estimation. The 
rotational invariance properties of both the transmit and receive arrays are investigated in [9], then the 
DOD and DOA are determined through two independent ID ESPRITs. However, an additional pairing 
operation is required. In [10], the relationship between two ID ESPRIT is investigated. In [11], the 
real-valued ESPRIT (unitary ESPRIT) is proposed to estimate DOD and DOA. It has lower 
computational complexity and slightly better angle estimation performance compared with 
ESPRIT [9,10]. A multi- singular value decomposition (multi-SVD) method is presented for DOD and 
DOA estimation in [12]. It provides better angle estimation than the traditional eigenvalue decomposition 
(EVD)/SVD method. The above schemes can only be used for angle estimation in the presence of 
spatial Gaussian white noise. In [13], an ESPRIT-based method for bistatic MIMO radar DOD and DOA 
estimation is proposed, which can eliminate spatial colored noise. However, it is only effective for three 
transmit antennas configuration. By dividing the transmit array into two subarrays, a combined ESPRIT 
and SVD of the cross-correlation matrix method (denoted as Chen's method) is presented in [14], which 
is effective for MIMO radar with three or more transmit antennas to eliminate the influence of spatial 
colored noise. 

However, in the subspace methods [13,14], the received signals are stacked into a special structure 
matrix, ignoring the multidimensional structure inherent in the received signals after matched filters. 
In this paper, a tensor-based frame is considered for the received signals, which exploits the 
multidimensional inherent structure and a novel tensor-based subspace for bistatic MIMO radar in 
the presence of spatial colored noise is proposed. Firstly, utilizing the multidimensional structure 
inherent in the received signals after matched filters, the received signals can be packed into a 
third-order measurement tensor. Then, the measurement tensor is divided into two sub-tensors, and a 
cross-covariance tensor is formulated to eliminate the spatial colored noise by exploiting the orthogonal 
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characteristic of matched filters. Finally, the higher-order singular value decomposition (HOSVD) 
technique is employed to formulate the signal subspace. The DOD and DOA are estimated through 
the ESPRIT algorithm, which are paired automatically. Theoretical analysis and simulation results 
validate that the proposed method suppresses spatial colored noise more efficiently and provides better 
angle estimation performance than Chen's method, the ESPRIT algorithm and the multi-SVD method, 
especially at the low signal-to-noise ratio (SNR) region. 

The rest of the paper is organized as follows. The tensor basics and signal model are presented 
in Section 2. A tensor-based subspace approach for angle estimation in the presence of spatial colored 
noise is proposed in Section 3. The computational complexity of the method is evaluated in Section 4. In 
Section 5, simulation results are provided to verify the performance of the proposed algorithm. Finally, 
Section 6 concludes this paper. 

AfotarfomScalars, column vectors, matrices and tensor are expressed by regular, bold lowercase, bold 
uppercase and bold calligraphic letters, respectively. [A]jj and [«4]j,j,jt stand for the (i, j) and k) 
element of a matrix, A, and a tensor, A. (-) H , (-) T , and (•)* denote the Hermitian transpose, 
transpose, inverse and complex conjugation without transposition, respectively. <g> and 0 denote the 
Kronecker operator and the Khatri-Rao product, respectively, diag(-) denotes the diagonalization 
operation, and arg(7) denotes the phase of 7. 

2. Tensor Basics and Signal Model 

2.1. Tensor Basics 

For the readers' convenience, several tensor operations are introduced firstly, which refer to [15,16]. 

Definition 1 (Matrix Unfolding): The three standard unfoldings of a third-order tensor, 
X e C IxJxK , denoted by [X] w G C IxJK , [X} {2) G C JxIK and [X] (3) G C KxIJ , can be expressed as 
[[<*1(i)k(fc-i)j+j = [X]i,j,k, [[X](2)]j,(i-i)K+k = [ x \i,j,k and [[X] {3) ] k:(j _ 1)I+i = [X] i:j:k , respectively. 

Definition 2 (Mode-n Tensor-Matrix Product):The mode-n product of X G C hxhx '" xlN by a 
matrix, A G C j " xl ", is denoted by y = X x n A, where y G c /lX/2X -' x/ "- lX J " x/ "+ lX - x/jv and 

In 
»n = l 

Definition 3 (The Properties of the Mode Product) :The properties of the mode product are shown 
as follows: 

X x n A x m B = X x m B x n A, m ^ n 

^x„Ax„B = ^x n (BA) (1) 

[X x 1 A 1 x 2 A 2 x--- x K A K ] (n) 
= A n .[A'] (n) .[A n+1 eg) A n+2 ... 0 A K ® Ax • • • ® A n _!] T (2) 

2.2. Bistatic MIMO Radar Signal Model 

Consider a narrowband bistatic MIMO radar system with M colocated antennas for the transmit array 
and N colocated antennas for the receive array, shown in Figure 1 . 
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Figure 1. Bistatic multiple-input multiple-output (MIMO) radar scenario. 




12 M 1 2 N 

Transmit array Receive array 



Both the transmit array and receive array are uniform linear arrays (UALs), and the inter-element 
spaces of the transmit and receive arrays are half-wavelength. At the transmit array, the transmit antennas 
emit the orthogonal waveforms S = [si, s 2 , • • • , sa/] t G C MxK , where K is the number of samples per 
pulse period. All the targets are modeled as a point-scatterer in the far-field, and it is assumed that there 
are P uncorrelated targets in the same range-bin of interest. {tp p }p =1 and {9 p }p =1 are the DOD and DOA 
with respect to the transmit and receive array normal, respectively. We consider a coherent processing 
interval (CPI) consisting of L pulses, then the baseband received signal for the /-th pulse period at the 
output of the receive array can be written as [3,14]: 

X, = BE { AS + W,, / = 1,2,...,L (3) 

where B = [b(0i), • • • , b(9 P )] G C NxP and A = [a(</?i), • • • , a.((pp)] G C MxP are the receive steering 
and the transmit steering matrix, respectively; h(0 p ) = [1, e jnsin9p , ■ ■ ■ , e MN-i)sme p jr e £Nxi and 
&(ip p ) = [1, e J7rsmVp , • • • , e i 7r ( M - 1 ) sm ^p] T g C A/Xl are the receive steering vector and transmit steering 
vector of the p-th target, respectively. £j = daig(c/) with c t = [f3 1 e j27T ^ dllTr , /3 P e j27T ^ dplTr ], f dP is 
the Doppler shift of the p-th target and T r is the pulse repetition interval and {(3 p }p =1 are the reflection 
coefficients. G C NxK is the noise matrix, and the columns of W; are independent and identical 
distribution complex Gaussian random vectors with zero mean and an unknown covariance matrix, 
Q. Unlike traditional phased-array radar, MIMO radar transmits mutually orthogonal waveforms, i.e., 
(l/fT)s m s^ = l,SjS? = 0(i,j = 1,2, ...M,i ^ j). Then, the received signals are matched by M 
transmitted waveforms, respectively. For the /-th pulse period, the output of the matched filter with the 
m-th transmitted waveform can be expressed as: 

Y ( , m = BD m cT + N iim) Z = 1,2,...,L (4) 

where Y^ m = (l/L)X z s^ G C^ 1 , D m = diag([a m (^i), a m (ip P )]), where a m ((f p ) is the m-th 
element of a(y? p ). N; m = (1/L)WjS^ G C Nxl is the noise vector after the matched filter with the 
m-th transmitted waveform, which is an independent, zero-mean complex Gaussian distribution with an 
unknown covariance matrix, Q, and satisfied with EfN^, N^] = 0(i,j = 1, 2, ...M, i ^ j). 
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3. Tensor-Based Subspace Approach for Angle Estimation 

According to Equation (4), the received signals are matched with all the transmitted waveforms. Then, 
we have: 

Y l = [Y M , Y hM ] = BEjA + N,, I = 1, 2, L (5) 

In the conventional subspace-based methods, the received signals in Equation (5) are packed into a 
special structure matrix as Y = [vec(Yi)), vec(Y 2 ), vec(Yx,)], which ignores the multidimensional 
structure inherent in the received signals. Based on the concept of Definition 1, it can be seen that 
the received signals for each pulse is a slice of a third-order tensor along the pulse direction. Therefore, 
by collecting L pulses, a third-order measurement tensor, y, is constructed, which is satisfied with 
[y]J = Y. According to E[N M , Ng] = 0(i,j = 1, 2, ...M, i ^ j), it is implied that the output of the 
colored noise has the orthogonal characteristic between different matched filters. In order to exploit this 
characteristic, the measurement tensor, y, is divided into two sub-tensors, which is shown as: 

yi = y*2 Fi (6a) 
y 2 = ^x 2 F 2 (6b) 

where Fi = [I Ml , 0 Mi x(m-M!)], F 2 = [0 M2X(M -m 2 ), Im 2 ] with M x + M 2 = M. According to 
Equation (6), we have: 

[tt]:,:,, = [Y u ,...,Y iMl ] = BE; A x + N] , I = 1,2,..,L (7a) 

[y 2 ],,i = [Y, Ml+ i, Y,, M ] = BS Z A 2 + N ; 2 , 1 — 1,2, L (7b) 

where Ax = F X A and A 2 = F 2 A. According to Equation (7), it is indicated that the measurement 
tensor data, and y 2 , are obtained from different matched filters, i.e, is the output of the first 
Mi matched filters and y 2 is the output of the residual M 2 = M — M x matched filters. Thus, we 
have E[(Ni) H ,Nf] = 0(1 = 1,2,...,L) owing to E[N M ,N«] = 0(i,j = 1,2,. ..M,i ± j). Then, a 
fourth-order cross-covariance tensor, H 2 i € C 7VxM2xArxMl , is formulated as: 

n 21 = ^y 2 . x (8) 

where [^i]n,,,ij = l/iEtiM^^i]^' = = l,-,M 2 .j = l,...,Mi. Since the 

spatial colored matrix, N/, and N 2 are also satisfied with E[(N i 1 ) H ,Nf] = 0(/ = 1,2,...,L), the 
influence of spatial colored noise is eliminated in Equation (8), i.e., the cross-covariance tensor, 1Z 2 i, 
is not affected by the additive spatial colored noise. According to Equation (8), the relationship between 
the cross-correlation matrix R 21 = ^^l^tVi]^) an d the cross-covariance tensor, H 2 \, is shown in 
Equation (9). 
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Then, the HOSVD [15,16] of the cross-covariance tensor, 1Z 2 i, yields: 

n 21 = s X1U1 x 2 u 2 X3U3 X4U4 (io) 

where 5 = ft 21 x x XJf x 2 x 3 x 4 G C WxM2XjVxMl denotes the core tensor [15] satisfying 
the property of all-orthogonality, while Ui,U 3 G C AfxJV , U 2 G C M2xMa and U 4 G C MlxMl are unitary 
matrices. Since 1Z 21 is a rank-P tensor, a cross-covariance subspace tensor, JF S , can be estimated by 
using the truncated HOSVD of 1Z 2 i, which can be written as: 

T s = <S S X1U1, x 2 U 2s x 3 U 3s x 4 U 4s (11) 

where U is (i = 1,2,3,4) contains the first P dominant singular vectors of Uj, S s = 1Z 2 \ Xi U^ s x 2 
x 3 U 3s x 4 U 4s as the reduced core tensor. Then substituting S s into Equation (11), we have: 

Fs =n 21 x 1 (U ls Uf s ) x 2 (U 2s U") x 3 (U 3s Vl) x 4 (U 4s Vl) (12) 

According to the relationship between the cross-correlation matrix and its 
corresponding cross-covariance tensor in Equation (9) and the Definition 3, a new 
cross -correlation matrix, R 2i , is reconstructed from J 7 ,, which can be expressed as: 

R 2 i = [(U ls UfJ ® (U 2s U 2 i s )]R 21 [(U 3s U 3 1 J ® (U^UfJ]* (13) 

In the subspace method [14], the signal subspace matrix, U s , is determined by using the truncated 
SVD of R 21 , i.e., R 21 xs U S A S V S . Inserting it into Equation (13), we have: 

R21 = {[(U ls Uf s ) ® (U 2s U 2 i s )]U s }A s {[(U 3s U 3 1 s ) ® (U 4s U 4 i J] T V s } H (14) 

According to Equation (14), using the truncated SVD of R 2i , the signal subspace, U s , can be written 

as: 

U s = [(U ls Ufj ® (U 2i U£)]U, (15) 

According to Equation (15), it is indicated that the signal subspace, U s , and XJ S span the same 
subspace. Hence, there exists a nonsingular matrix, T G C PxP , satisfied with tj s = (A 0 B)T. 
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After obtaining the signal subspace, U s , the ESPRIT algorithm [7,8] is applied to estimate the DOD 
and DOA. 

In order to estimate both the DOD and DOA, the signal subspace, U s , is divided into four submatrices: 
U sl = IW,, U s2 = r 2 tJ s , U s3 = r 3 tJ s and U s4 = r 4 tJ s , where T l = J^ 2 " 1 * <g> 1^, 

1*2 = J(^ 2 ^ ® Iat, r 3 = 1m 2 <S> J(^) X \ r 4 = Im 2 <S> 1 \ = [I& 0fc X i], j[ 2 j = [Ofcxi In 
doing so, we have: 

T 2 U S = TiU^ (16a) 

r 4 u s = r 3 u s * r (i6b) 

where * t = T ^/T, * r = T ^/T, $ t = diag([e^ sin ^, • • • , e jnsin ^ p ]) and $ r = diag([e^ sin91 , 
• • • , eJ'Tsinep]^ con tain the desired DOD and DOA information. Equation (16) can be solved by least 
squares (LS), then, and \I> r are obtained. Let <fr t and T be the eigenvalue matrix and eigenvector 
matrix of ty t - Then, the DOD of the p-th target is derived as: 

ip p = arcsin(arg(7^)/7r), p=l,2,...,P (17) 

where 7^ is the p-th diagonal element of <f> t . Note that \I> t and \I> r have the same eigenvector matrix, the 
diagonal matrix, $ r , can be determined as $ r = T _1 Vl/ r T. Then, the diagonal elements of $ t and <3? r in 
the same position correspond to the same target, i.e., the DOD and DOA are paired automatically. The 
DOA of the p-th target is derived as: 

6 P = arcsin(arg(7^)/vr), p = l,2,...,P (18) 

where 7^ is the p-th diagonal element of <fr r . 

4. Computational Complexity Analysis and Remark 

In order to analyze the computational complexity of the proposed method, it is necessary to know the 
complexity of the SVD algorithm. There are a lot of methods to compute SVD, and the computational 
complexities of them are different. In [17], it has been pointed out that orthogonal iteration is an efficient 
solution for SVD algorithm. The computational complexity of an M x N matrix truncated to rank r is 
0(k r MNr) by using this orthogonal iteration, where k r is a constant that depends on the design of the 
algorithm. The main computational burden of the proposed method, multi-SVD method, Chen's method 
and ESPRIT algorithm is the estimation of the signal subspace. In order to estimate the signal subspace, 
tj s , the proposed method needs to calculate the truncated HOSVD of 1Z 2 i and the truncated SVD of 
R21. The truncated HOSVD of 1Z 2 i is equivalent to the truncated SVD of all its matrix unfolding, which 
needs 0(Ak r MiM 2 N 2 P). Additionally, the computational complexity of the truncated SVD of R 24 is 
0(k r M l M 2 N 2 P). The total computational complexity of the proposed method is 0(5k r M 1 M 2 N 2 P). 
In the multi-SVD method, the signal subspace is estimated by the truncated HOSVD of the third tensor, 
y. Thus, the computational complexity is 0(3k r MNLP). Chen's method only needs the truncated 
SVD of R 2 i to estimate the signal subspace, which needs 0(£; r MiM 2 iV 2 P). The ESPRIT algorithm 
uses the truncated SVD of the covariance matrix R = (1/L)YY H to estimate the signal subspace, which 
needs k r M 2 P. According to the above analysis, the computational complexity of the proposed method 
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is similar to the multi-SVD method, but higher than both Chen's method and the ESPRIT algorithm. 
However, the proposed method provides better angle estimation performance than all the aforementioned 
methods, which is demonstrated in the next section. 

Remark 1: According to Equation (11), in order to obtain the signal subspace, XJ is , of 
Uj (i = 1, 2, 3, 4), the necessary conditions for M lt M 2 , N and L are that: M 1 > P, M 2 > P, N > P, 
L > P. It is indicated that the maximum number of targets can be identified by the proposed method 
is min[Mi, M 2 , N, L\. Thus, the number of targets that can be identified by the proposed method is 
smaller than ESPRIT algorithm. However, the proposed method performs well and provides better angle 
estimation performance in the presence of spatial colored noise, while the ESPRIT algorithm has marked 
performance degradation, especially in the low SNR region. 

5. Simulation Results 

In this section, some simulations are presented to evaluate the angle estimation performance of the 
proposed method in the presence of spatial colored noise. The multi-SVD algorithm [12], ESPRIT 
algorithm [10] and the method in [14] (denoted as Chen's method) are used to compare with the 
proposed method. We consider a MIMO radar system with M colocated antennas and N colocated 
antennas for the transmit and receive array, respectively. Both of the transmit array and receive array 
are half- wavelength spaced ULAs. M transmit antennas transmit M orthogonal waveforms, and the 
m-th transmitted waveform is the m-th row of S G C KxK , where S = (1 + j)/y/2H.K, an d is 
the K x K Hadamard matrix. The number of samples per pulse period is K = 256, and the pulse 
repetition interval is T r = 5 us. There exists three uncorrected targets located at (<p 1: 9i) = (30°, —30°), 
(<P2,9 2 ) = (-40°, 10°) and (<p 3 ,0 3 ) = (10°, 10°), and the reflection coefficients of the targets are 
{P P Y P =i = 1- The Doppler shifts are {f dp } 3 p=1 = {300, 400, 500}Hz. The spatial colored noise 
is modeled as a spatial complex autoregressive (AR) model of second-order with the coefficients 
z — [1, —1, 0.2] [12], and the root mean square error (RMSE) of the angle estimation is defined as: 



the i-th Monte Carlo trial. 

Figure 2 shows the RMSE of angle estimation of different methods versus SNR, where M = N = 12, 
Mi = 3, L = 100 and Q = 200. It is shown in Figure 2 that the ESPRIT and multi-SVD algorithm 
provide worse angle estimation performance than Chen's method and the proposed method, especially 
at the low SNR region. This is because the ESPRIT and multi-SVD method cannot eliminate the 
influence of spatial colored noise. It also can be observed that Chen's method provides better angle 
estimation performance than the ESPRIT algorithm, which is consistent with [14]. Owing to taking the 
multidimensional structure into account and using the cross-covariance tensor technique, the proposed 
method can eliminate the spatial colored noise more efficiently. Thus, the proposed method outperforms 
all aforementioned methods, especially at the low SNR region. 




(19) 



where Q is the number of Monte Carlo trials and (p { and Q { are the estimation of DOD ipi and DOA 9i of 
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Figure 3 shows the probability of the successful detection of the different methods versus SNR, where 
M = N = 12, Mi = 3, L = 100 and Q = 200. Successful detection requires that the absolute 
error of both DOD and DOA for all three targets are within 0.5°. It can be seen from Figure 3 that 
all the methods exhibit a 100% successful detection at high SNR values. As the SNR decreases, the 
probability of successful detection decreases for each method at a certain point, which is known as the 
SNR threshold. It also can be seen that the proposed method has a lower SNR threshold than Chen's 
method, ESPRIT and the multi-SVD method, owing to the super capability of eliminating the spatial 
colored noise. 

Figure 2. Root mean square error (RMSE) versus SNR for P = 3 targets. 




-25 -20 -15 -10 -5 0 5 10 15 

SNR (dB) 



Figure 3. Probability of successful detection versus SNR. 




SNR (dB) 

Figure 4 shows the RMSE of angle estimation of different methods versus the number of pulses, where 
M = N = 12, Mi = 3, SNR= -5dB and Q = 200. It can be seen in Figure 4 that the angle estimation 
performance of all methods can be improved with pulse increases, and the proposed method provides 
better angle estimation performance than Chen's method, ESPRIT and the multi-SVD algorithm. 
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Figure 4. RMSE versus pulses for P = 3 targets. 
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6. Conclusions 

In this paper, a tensor-based subspace approach is presented to DOD and DOA estimation for bistatic 
multiple-input multiple-output (MIMO) radar in the presence of spatial colored noise. The proposed 
method exploits the the multidimensional structure inherent in the received signals to construct a 
third-order measurement tensor. Then, two sub-tensors are obtained from the measurement tensor, which 
can be used to formulate a cross-covariance tensor for eliminating the influence of spatial colored noise. 
Finally the DOD and DOA can be estimated in conjunction with the ESPRIT method. The proposed 
method has better angle estimation performance than Chen's method, ESPRIT and the multi-SVD 
method, especially at the low SNR region. Several simulation results have verified the performance 
of the proposed method. 
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